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Predicting the Merger Fraction of Lyman alpha Emitters 
from Redshift z ^ 3 to z ^ 7 
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ABSTRACT 

Rapid mass assembly, likely from mergers or smooth accretion, has been predicted to 
play a vital role in star-formation in high-redshift Lya emitters. Here we predict 
the major merger, minor merger, and smooth accreting Lya emitter fraction from 
« 3 to z w 7 using a large dark matter simulation, and a simple physical model 
that is successful in reproducing many observations over this large redshift range. The 
central tenet of this model, different from many of the earlier models, is that the star- 
formation in Lya emitters is proportional to the mass accretion rate rather than the 
total halo mass. We find that at z w 3, nearly 35% of the Lya emitters accrete their 
mass through major (3:1) mergers, and this fraction increases to about 50% at z w 7. 
This imply that the star-formation in a large fraction of high-redshift Lya emitters 
is driven by mergers. While there is discrepancy between the model predictions and 
observed merger fractions, some of this difference {^ 15%) can be attributed to the 
mass-ratio used to define a merger in the simulation. We predict that future, deeper 
observations which use a 3:1 definition of major mergers will find ^ 30% major merger 
fraction of Lya emitters at redshifts ^ 3. 

Key words: cosmology - theory: dark matter, galaxies-high-redshift: interactions: 
halos, methods: numerical 



1 INTRODUCTION 

While several theoretical models of Lya emitters 
have been developed (e.g. iBarton et al.l |2004|: iDave et al 



20061 ; lTasitsiomill2006l:IShimizu et al.ll2007l: iNasamine et al 



2008"; 'Lc Dclliou et al.' '20061; iKoba vashi et al.n2007l . l2009l : 
Daval ot al. 2008; Samui et al 2009) we still lack a complete 
understanding of how mass assembly and star-formation oc- 
curs in these galaxies. It is likely that mergers are the dom- 
inant mode of mass accretion resulting in star-formation 
in at least so me fraction (^ 20%) of Lya emitters (e 
iPirzkal et all boOT; Bond ct al. 2009; Taniguchi et al.ll200! 

Recently, iTilvi et al.l ((20091) developed a physical model 
of Lya emitters which is successful in explaining many 
observable including number density, stellar mass, star- 
formation rate, and clustering properties of Lya emitters 
from 2 « 3 to 2 « 7. This model d iffers fundamentally from 
many of the earlier models (e. 



Diikstra et aTl l2007l: 



Fernandez fc Komats' 



i Mao et al 
|2Q0^ 



Haiman fc Spaans 
l2007l ; IStark et al.l 



199S 



2007 



in that the star-formation of 
Lya emitters, and hence their Lya luminosity is propor- 
tional to the mass accretion rate rather than the total halo 
mass. 
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In this paper, we combine the above model with a large 
dark matter cosmological simulation to predict the major 
merger, minor merger, and smooth accreting Lya emitter 
fraction from z ~ 3 to ^ « 7. We also carefully asses the 
uncertainties in our model predictions. We note that cur- 
rently there are no theoretical prediction of merger fraction 
of Lya emitters. 

On the observational front, only recently it has been 
possible to estimate the m erger fraction of Lya emitters. At 
z ~ 0.3. ICowie et al.l JJOIQ) found that ^ 30% of Lya emit- 
ters are either irregulars or have disturbed morphologies in- 
dicative of mergers. At higher redshifts, z > 3, the observed 
merger fraction varies f rom » 20% to » 45% (jPirzkal ct ajj 
I2OO7I : iBond et al.|[2009l; iTaniguchi et al.|[2009l ) . 

Recently, I Cooke et al.l l|2010l ) found that all Lyman- 
break galaxies (LBGs) that have close companions (^ I5h~^ 
Mpc) exhibit strong Lya emission lines. In the spectro- 
scopically confirmed close pairs, the spectra/imaging show 
double Lya emission peak/double morpholog y confirming; 
that these LBGs at z « 3 are indeed mergers. ICooke et al 



(l20ld ) also studied a sample of spectrosco pically confirmed 
LBGs at z « 3 from IShaplev et al.l (|2006[ ). and found that 
about 33% of LBGs with strong Lya emission have double 
Lya emission peaks, strengthening the above conclusion. 
The outline of this paper is as follows. In §2 we briefly 
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descr ibe the physical model of Lya emitters IjTilvi et alj 
l2009r ). our simulations, and how we construct the merger 
trees. In §3, we estimate the merger fraction of Lya emit- 
ters, discuss the effect of mass resolution in the simulation 
on our results, and compare our model predictions with the 
observations. In §4 we quantify the uncertainties associated 
with the predictions, and we summarize our results in §5. 



2 METHODS 

2.1 Modeling Lya Emitters 

Our method is based on carrying out dark matter (DM) 
only simulations, and then populating the halos from these 
simulations with Lya emitters accordin g to a simple physi- 
cal model presented in lTilvi et al.l (|2009l ) . This model uses a 
single parameter to successfully reproducing many observed 
properties of Lya emitters, including their luminosity func- 
tions, stellar masses, stellar ages, star formation rates, and 
clustering from z — 3.1 to 2 = 6.6. The central tenet of this 
model is that the Lya luminosity is proportional to the 
star formation rate i.e. 



JLya 



= 1 X 10* 



SFR 

Meyr- 



erg s 



(1) 



which in turn is proportional to the mass accretion rate i.e. 
SFR =f* X Mb. Here, f* and Adt are the star formation effi- 
ciency (i.e. converting baryonic mass to stars), and the bary- 
onic mass accretion rate. As we use DM only simulations, 
the baryonic mass accretion rate is obtained by converting 
DM mass accretion assuming universal ratio of baryonic to 
DM densities . For more detail s about this model, we refer 
the reader to lTilvi et all (|2009D . 



2.2 Simulation and Halo Catalogs 

In order to implement this method, we carried out a large 
N-body DM cosmological simulation (hereafter Gadget- 
1024), with the GADGET2 (Springcl 2005) code. We gener- 
ated the initial conditions for the simul ation using second- 
order Lagrangian Perturbatio n Theory (|Crocce et al.l 120061 : 
iThacker fc Couchman I [200a). In this simulation we use 
1024^ DM particles, each particle with a mass Mp r; 2.7x 10^ 
Mq/i"^ evolved in a comoving volume of (102 Mpc)^ . Usin g 
a Friends-of-Friends (FOF) halo finder l|Davis et al.lll985h . 
we identify DM halos that contain 100 or more DM particles. 
This corresponds to a minimum halo mass Mhaio ~ 2.7 x 10^ 
Mq/i-i. 

We then generate catalogs, for redshifts from z = 10 
to 2; = 3, which contain positions of halos, their DM mass, 
and unique IDs of each individual particle that belongs to a 
given halo. These unique particle IDs are later used to track 
halos between two epochs. Throughout this chapter we as- 
sumed a flat ACDM cosmology with parameters f2m=0.233, 
nA=0.721, ^^6=0.0462, /i=0.71, a8=0.817 where f7„, Qa, 
fib, h, and ag correspond, respectively, to the matter den- 
sity, dark energy density, and baryonic density in units of 
the critical density, the Hubble parameter in units of 100 
km s~^ Mpc"'^, and the RMS density fluctu ations on the 
8 Mp c h~^ scale, in ag reement with WMAP (jSpergel et al.l 
120071 ) flve year results (JHinshaw et al.ll2009l '). 
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Figure 1. Merger schematic showing classification of mergers 
based on the progenitor mass ratio. 



2.3 Lya Emitter Catalogs 

To construct a simulated population of Lya emitters, we 
first calculate the amount of D M mass that is accreted by 
each halo in time tLya=SO Myr JTilvi et al.ll2009l ). This time 
interval was chosen so as to match with the average stel- 
lar population of Lya emitters. This mass accretion rate 
is then converted to the Lya luminosity as described in 
Section 2.1. While, in general, we expect every halo to ac- 
crete more mass with time, we find that some halos lose 
mass (negative accretion) which is due to the limitation of 
halo-finding technique. We correct for this simulation noise 
by counting the number of halos that have negative mass 
accretion in each mass bin and subtracting this from corre- 
spond ing counts in the positive mass bins (See iTilvi et al.l 
(200!?) for more details). With this procedure we now have 
Lya emitter catalogs from z = 6.6 to 2 = 3. 

It is possible that the mass resolution in our simu- 
lation might affect our results. In order to carefully as- 
sess this effect, we also us ed the Millenium-II simulation 
(|Bovlan-Kolchin et al.ll20 0g) which has a larger volume, and 
a higher halo mass resolution. In particular, it follows 2160'' 
particles with each particle having a DM mass 6.885 x 10® 
M0/i~^ in a simulation volume of (141 Mpc)''. Further- 
more, in this simulation, each halo is identified by an unique 
ID allowing us to follow the same procedure as applied to 
our Gadget- 1024 simulation. We use Millenium-II simula- 
tion only to investigate the simulation mass resolution effect 
since the time step between each output is about ten times 
larger, making it unsuitable for our final comparisons with 
observations. 



2.4 Merger Tree 

As described in Section 2.1, each Lya emitter is assigned 
a luminosity based on its mass accretion rate. This mass 
accretion can occur when either two or more progenitors 
(halos at time step S-1) merge into a single descendent (halo 
at a later time step S), or when a single progenitor evolves 
into a single descendent (see Figure 1). 

We then classify each merger into one of three cate- 
gories: major merger, minor merger, or smooth accretion, 
based on the progenitor mass ratio. In particular, we define 
each of these cases as: 

(1) major merger: 0.3 ^ m2/mi ^ 1, 

(2) minor merger: 0.1 ^ m2/mi < 0.3, and 

(3) smooth accretion: m-ilni\ < 0.1, 
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Figure 2. The fraction of Lyo emitters that have formed from 
major mergers, minor mergers, and smooth accretion at 2 = 3.1. 
The major and minor merger fractions are shown in soHd, and 
dotted lines respectively, while the long dashed line shows the 
fraction of Lya emitters that accrete mass through smooth 
accretion. The shaded area are the poisson errors on the major 
merger fraction. About 35% of the LyQ emitters at z i5i 3 accrete 
their mass through major mergers. The steep decline in the major 
merger fraction at Log Lya < 42 erg s~^ is not real but results 
from limited simulation mass resolution (see Section 3). 



where m\ and m2 are the masses of the most massive, and 
second-most massive progenitors associated with a single de- 
scendant (see Figure 1). In the case where a descendent halo 
has a single progenitor halo, by defauh, this halo is assigned 
a smooth accretion mode. Using the above criteria we now 
have a merger tree, i.e. each descendent halo, and hence 
each Lya emitter, has been associated with the mode of 
mass accretion i.e. major merger, minor merger or smooth 
accretion. 



3 MODEL RESULTS AND OBSERVATIONAL 
COMPARISONS 

3.1 Merger Fraction at 2: « 3 

Using this method we are able to calculate the fraction of 
Lya emitters that have undergone major mergers, minor 
mergers, and smooth accretion. Here, we remind the reader 
that each descendent halo hosts one Lya emitter. In effect, 
every descendent halo is an Lya emitter with its luminos- 
ity proportional to the mass accretion rate. We define the 
merger fraction as the ratio of number of Lya emitters 
formed from mergers in a given luminosity bin to the total 
number of Lya emitters in that bin. Thus, 

Nm{L) 



fiL) = 



Nt{L) 



(2) 



where Nm{L) is the number of Lya emitters formed from 
mergers, and Nt{L) is the total number of Lya emitters 
in that luminosity bin. 

Figure 2 shows the merger fractions of Lya emitters 
that have formed from major mergers, minor mergers, and 
smooth accretion, at z — 3.1. The major merger, and minor 
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Figure 3. Comparing the merger fractions between Gadgetl024 
and the Millenium-II simulation. The vertical dotted lines repre- 
sent the luminosity range below which our results are not reliable. 
In the Millenium-II simulation, due to higher mass resolution, the 
vertical dotted line shift towards the lower luminosity. Also seen is 
the decreased noise at the brighter luminosities in the Millenium- 
II simulation due to its larger simulation volume. 



merger fractions are shown in solid and dotted lines respec- 
tively, while the long dashed line indicates the fraction of 
Lya emitters formed from smooth accretion. The shaded 
region shows the Poisson errors on the major merger frac- 
tion. 

From Figure 2 it is evident that at z = 3.1, about 35% 
of Lya emitters are formed from major mergers. In other 
words the dominant mode of mass accretion, and hence star- 
formation in nearly 35% of Log(Lya) > 42 erg s~^ emit- 
ters, occurs through major mergers, while remaining 65% 
of Lya emitters accrete their mass through minor mergers 
and smooth accretion combined. 

Note that the Poisson errors dominate the bright end of 
the Lya luminosity. Our results might not be fully reliable 
below 6 X 10*^ erg s~^ where we see a steep decline in the ma- 
jor merger fraction. This effect might be partly due to lower 
halo mass resolution. To understand and quantify this effect 
of lower halo mass resolution, we now compare our results 
from figure 2 with the Millenium-II simulation which has 
about 2.6 times more volume, and nearly forty times better 
mass resolution. In addition, we also investigate the effect 
of stellar ages ( tLya) of Lya emitters on our predictions. 



3.2 Effect of Halo Mass Resolution and tLya 

Using the same procedure as described in Section 2, we 
constructed a Lya emitter catalog from the Millenium- 
II simulation, at a common redshift, z = 6.8 between the 
two simulations. In Figure 3 we compare the results from 
the two simulations. The top panel are the results from our 
Gadget- 1024 simulation while the bottom panel shows the 
merger fraction obtained from Millenium-II simulation. The 
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II,. ~ 30 Myr 
tL,a ~ 60 Myr 
II,. ~ 100 Myr 




The major merger fraction (left panel), and smooth accre- 
tion (right panel) show a mild evolution with redshift. On 
the other hand, the minor merger fraction is nearly constant. 



42.5 43.0 

Lya Luminosity (erg s~^) 



Figure 4. Model prediction of major merger fraction for different 
stellar ages of Lya emitters at 2 ^ 7. The smaller t^^Q implies 
younger stellar population in Lya emitters. 



vertical dotted lines show the luminosity below which the 
apparent decline in the major merger fraction is caused due 
to lower halo mass resolution. This luminosity range shift 
towards lower Lya luminosity in the Millennium-II simu- 
lation (lower panel) which has a better mass resolution. Also, 
at the bright end of the luminosity axis, the fluctuation in 
the merger fraction is lesser in the Millenium-II simulation 
due to the larger volume. Thus, better mass resolution is 
needed to estimate the merger fractions reliably at lower lu- 
minosities, and a larger simulation volume to minimize the 
noise at the brighter luminosities. Hereafter, in our analy- 
sis, we consider only those Lya emitters that have Lya 
luminosity > 6 x lO^^erg s~^. 

We now study the effect of t^yaon our model predicted 
results. For our model Lya emitters the luminosity is pro- 
portional to the mass accretion rate (AAfj, / thy a) where 
we assumed a fixed t_Lya=30 Myr which corresponds to the 
average stellar ages of Lya emitters. While this is true 
for average population, from observations it is inferred that 
this stellar age ranges f rom 1 Myr to even hundred Myr (e.g. 
iFinkelstein et al.ll2007l '). 

Figure 4 shows the major merger fractions of our model 
Lya emitters a,t z = 6.6 for tLya^SO, 60, and 100 Myr. 
While the stellar ages have changed by a factor of 3, we 
find no or little change in major merger fraction. Thus, our 
model predictions are independent of the stellar ages of Lya 
emitters. 



3.3 Redshift Evolution of Merger Fraction 

In preceding sections we have shown that the predicted 
merger fractions are affected by the mass resolution, and 
that there is no dependence of tiy^on the predicted merger 
fractions. We now investigate the redshift evolution of 
merger fraction of Lya emitters. 

Figure 5 shows the redshift evolution of major merger 
fraction (left panel), minor merger fraction (middle panel), 
and smooth accretion (right panel) from z = ?,.l to z = 6.6. 



At z = 3.1 the average major merger fraction /„ 



; 35%. 



Similarly, the average fraction of minor mergers (middle 
panel) and smooth accretion (right panel) are f minor ~ 25%, 
and fsmooth ~ 40% respectively. At z = 6.6, the major 
merger fraction is higher with f major ~ 50%. Thus, at higher 
redshifts, the mass accretion and hence the star-formation 
in a significant fraction of the Lya emitters occurs through 
major and minor mergers. 

3.4 Comparison with the Observations 

In this section we compare only major merger fractions with 
the observations since it is extremely difficult to identify mi- 
nor mergers from observations especially at higher redshifts. 

Only recently, has it been possible to estimate the ob- 
served merger fraction of Lya emitte r s at h igh-redshifts. 
For example, at z — 5.7 lTaniguchi et al.l (|2009l ) have studied 
morphological properties of Lya emitters in the COSMOS 
field using HST/ACS. Only two out of about 50 Lya emit- 
ters show clear signs of extended features, an indication of 
either interacting or merging galaxies. However, there is a 
large uncertainty in estimating merger fraction from these 
observations due to their shallow survey depth. 

In Figure 6, we compare this result with our model pre- 
dictions. The downward arrow indicates the upper limit on 
the merger fraction at z = 5.7. At this redshift, our model 
predicted major merger fraction is about 50%. While com- 
paring our model predictions with the observations we have 
accounted for the observed limiting Lya luminosity at cor- 
responding redshifts. 

At slightly lower redshift z ^ 5, (jPirzkal et all 120071 ') 
studied morphologies of nine Lya emitters in HUDF. 
To quantify the morphologies, th ey used concentration (C) , 
and asymmetry{A) parameters (|Conselice et al.ll200d ) and 
found that nearly 44% have clumpy or complex structures. 
Visually, about 33% sources look morphologically disturbed 
or as ongoing mergers, while nearly 10% of the sources can 
not be reliably identified as mergers. This observed merger 
frac tion is nearl y same as our model prediction. 

iBond et al.l (|2009l ) studied morphological properties of 
about 120 z ~ 3.1 Lya emitting galaxies in the rest-frame 
ultra-violet band. They found that at least 17% of the total 
Lya emitters contain multiple components which might in- 
dicate that these are either individual star-forming regions 
within a single galaxy, a merged system or ongoing merg- 
ers. Since it is very difficult, due to their compact sizes (e.g. 
Malhotra et al 2011, Bond et al 2011), to definitely con- 
clude whether the multi-component Lya emitters are the 
remnants of mergers or if these are individual star-forming 
regions in a single system, we have shown the merger fraction 
by upward and downward arrows indicating uncertainties in 
both directions. Their classification is based on counting the 
number of components in a fixed aperture. At z — 3.1, our 
model predicted major merger fraction is about 35%, much 
higher than the observed fraction. However, in Section 4 we 
show that, some of this difference between model predictions 
and observations can be attributed to the definition of major 
merger mass ratio. 

At lower redshift, z ^ 0.3, ICowie et all (|2010l ') stud- 
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Figure 5. Redshift evolution of major merger fraction (left panel), minor merger fraction (middle), and smooth accretion (right panel) 
from z = 3.1 to ^ = 6.6. We see a mild evolution of major merger fraction (left), and in smooth accretion (right). On the other hand the 
minor merger fraction remains constant over this redshift range. 



led morphologies of Lya emitters in the GROTHOO and 
SIRTFFLOO fields. They found that > 30% of the Lya 
emitters show signs of ongoing mergers. Based on the above 
comparisons, there is some discrepancy between our model 
prediction and observations. In the following section, we in- 
vestigate the uncertainty due to the mass ratio used to define 
a merger. 



4 UNCERTAINTIES FROM MODEL 

PREDICTIONS AND OBSERVATIONS 

4.1 Dependence of Merger Fraction on the 
Progenitor Mass Ratio Definition 

We now vary the progenitor mass ratio criteria from 1:3 to 
1:2, i.e. major mergers are now defined as those for which 
0.5 ^ 77i2/mi s; 1, while for 0.1 ^ m2/m\ < 0.5, they are 
classified as minor mergers. All other Lya emitters with 
their progenitor halo mass ratio m2/m\ < 0.1 are defined 
as smooth accreting. In figure 7 (see also Figure 2) we show 
the predicted merger fraction dependence on the progenitor 
mass ratio definition. 

By changing the progenitor mass ratio from 1: 3 to 1: 
2, the major (minor) merger fraction drops (increases) by 
about 15%. Thus, it is clear that the predicted merger frac- 
tion of Lya emitters depends on the the progenitor mass 
ratio definition, and one needs to be careful when compar- 
ing model predicted merger fractions with the observations 
since the merger defining mass ratio infiuences the merger 
fractions. 



4.2 Merger Observability Timescale 

In addition to the above uncertainty, the observed merger 
fraction depends on the timescale, tabs (observability 
timescale) during which a merger can be identified as a 
merger. For example, if tabs is shorter than tLya, the time 
during which a galaxy is observed as an Lya emitter, 
then the observed merger fraction will be underestimated. 
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Figure 6. Comparison of our model predicted major merger frac- 
tion with the observations. Filled circles are our model prediction 
while other symb ols represent observ ations at z = 5.7 (Taniguch i 
et al 2009), z ^ 5 llPirzkal et al.l2o"o7l l. z = 3.1 l|Bond et al.l2009h . 
and 2 Ri 0.3 jCowie et al.ll201CI) . Error bars on model predictions 
indicate Poisson errors. 



According to our model (see equation 1), merging of two 
galaxies results in an Lya emitter with its Lya lumi- 
nosity proportional to the mass accretion rate. If, say two 
galaxies merge on a dynamical timescale (fdyn), then it will 
be observed as a merger for a period tdyn=tobs. 

Using tdyn = •^37r/(16Gp(z)), where p{z) is the average 
overdensity of the collapsing object, we find that at z «3, 
tdyn ~ 300 Myr, much larger than tLi;a=30 Myr. Thus, we 
expect to miss no or little merger fraction of Lya emitters 
at z ^ 3, due to the observability timescale. The dynamical 
time, however is nearly half with tdyn ~ 130 Myr at a « 6. 

As can be seen from the above discussion that while the 
predicted merger fraction is not affected by the merger ob- 
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assembly and hence the star-formation in majority of the 
Lya emitters. 

• We also compared our model predicted major merger 
fraction with the observations and found that our model 
over-predicts the major merger fraction. This discrepancy 
can be resolved if we take into account the uncertainties 
in defining the merger classification mass ratio. By chang- 
ing the major merger mass ratio from 3:1 to 2:1, we found 
that the predicted major merger fraction decreases by about 
15%. We predict that future, deeper imaging surveys and us- 
ing spectroscopic m ethods such as the one demonstrated in 
ICooke et al.l ([20101), should find ^ 30% merger fraction of 
Lya emitters at high-redshifts. 

This work was supported in part by the grant HST- 
0808165. All simulations were performed on the saguaro 
cluster operated by the Fulton School of Engineering at Ari- 
zona State University. 



Figure 7. Progenitor mass ratio dependence on the predicted 
merger fraction of Lya emitters at 2: = 3.1. Here, the major 
merger is defined as 0.5 ^ m^/mi ^ 1, while for the minor merger 
0.1 Sj m2/mi < 0.5. All Lya emitters with its progenitor mass 
ratio ra2lra\ < 0.1 are defined as smooth accreting. 



servability timescale, there is some uncertainty in the pre- 
dicted merger fraction due to the defining mass ratio. In 
addition, a large uncertainty is due to the shallower survey 
depths especially at higher-redshifts. 



5 SUMMARY AND CONCLUSIONS 



Motivated by our earlier work (jTilvi et al.l 12009 ) based on 
a simple physical model of Lya emitters, combined with 
a large dark matter cosmological simulation, we have pre- 
dicted the major merger, minor merger, and smooth accret- 
ing Lya e mitter fraction from z ~ 3 to 2: ~ 7. The model 
presented in iTilvi et al.l ((20091), different from many of the 
earlier models in that the star-formation rate is proportional 
to the mass accretion rate rather than the total halo mass, 
has been successful in reproducing many observed physical 
properties including the luminosity functions, stellar ages, 
star formation rates, stellar masses, and clustering proper- 
ties of Lya emitters from 2; « 3 to z ~ 7. 

We carefully constructed the merger tree from z « 7 to 
z « 3, accounting for the uncertainties in the halo finder. In 
this merger tree, each descendent halo was associated with 
its progenitor(s), and based on the progenitor mass ratio, 
each merger was classified as either major, minor merger or 
smooth accretion. We also carefully assessed how our pre- 
dicted results might be affected by several parameters in- 
cluding halo mass resolution, defining merger mass ratio, 
stellar ages of Lya emitters, and the merger observability 
timescale. We summarize our key results below: 

• At z « 3, about 35% of Lya emitters accrete their 
mass through major mergers, and this fraction increases to 
about 50% at redshift 2 « 7. On the other hand, at this 
redshift, the minor mergers and smooth accretion contribute 
equally in the remaining 50% of the Lya emitters. Thus, 
at higher redshifts, mergers play an important role in mass 
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